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PREFACE 


This document is submitted by the McDonnell Douglas Astronautics Company 
to the National Aeronautics and Space Administration and was prepared 
under Contract NAS8-33191, “Statistical Energy Analysis of Complex Struc 
tures (Phase 2),” The study was directed by R. W. Trudell. W. Clever 
of the Vibration Analysis Branch, Systems Dynamics Laboratory, Marshall 
Space Flight Center, administered and directed the contract. 
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SYMBOLS 


area 

surface area of element exposed to acoustic excitation 

speed of sound in air 

group velocity = 1.07[ujCj^t] 

longitudinal wave velocity 

Eh^ 

plate bending rigidity = t 

12(l-v ) 

modulus of elasticity, total energy of element 

joint length 

number of modes 

power input 

accel erati on 

frequency 

ring frequency 

gravitational acceleration 

plate thickness 

radius of gyration of plate 

length 

mass 

modal density 

pressure 

radius 


thickness 


V 



velocity 
weight density 

SEA response equation coefficient 
A increment 

n element loss factor 

Y Poisson's ratio 

TT 3.14159 

p mass density 

a radiation efficiency 

T joint loss coupling parameter 

4,.. power transfer coefficient for coupling between modes in elements 

i and j 

a> frequency in radians/second, normally center frequency of a frequency 

band 



Element 

System 



a set of modes modeled as one unit of a system 

the total structure and associated energy sources under consideration 
indicates averaging over both time and space 
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Section 1 
INTRODUCTION 

Statistical Energy Analysis (SEA) is a method of estimating the structural 
vibration properties of complex systems in high frequency environments. The 
method considers the average distribution and transfer of energy among the 
modes of a vibrating system and assumes that all of the vibratory energy of 
the systan is contained in these modes. The analysis method generally 
involves averaging structural responses over time and space. 

The response predictions are made by modeling the structure into relatively 
gross elements, deriving power flow equations between these elements and the 
environment, and solving the resulting system of equations for element response 
1 evels. 

The general assumptions which SEA is based on are: 

1) The total vibrating system can be partitioned into SEA elements 
(with suitable boundary conditions) whose modes approxima^^ft 
modes of the original vibrating systan. 

2) The modes of the elements of a system contain all of the vibratory 
energy of the system. 

3) Only modes occurring within the same frequency band are coupled. 

4) The energy in one frequency band of a system element is equally 
distributed among the modes of that element occurring in the 
frequency band. 

5) For two coupled elements, all of the modes occurring in one of the 
elements in one frequency band are equally coupled to each mode 
occurring in the same frequency band in the other element. 

A more thorough treatment of SEA is contained in the references. 
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Section 2 

VIBRATION PREDICTION FOR PAYLOAD/SUPPORT STRUCTURE 
WITH ACOUSTIC EXCITATION 

The structure analyzed was the Materials Experiment Assembly (MEA) (Fig- 
ure 1), a portion of the OSTA-2 payload for the Space Transportation 
System, SEA techniques were used to model the structure and predict 
structural element responses to given acoustic excitation. The predicted 
responses were compared with the physical test results obtained at MSFC 
for evaluation of SEA accuracy. 



Figure 1. OSTA-2 Aft MEA MPE Support Structure 
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2.1 MODELING 


The structure was represented by six SEA elements chosen, to the maximum 
practical extent, to correspond wich MEA acoustic test accelerometer 
locations. This was done to assure valid comparisons between the response 
predictions and the corresponding acoustic test data. A schematic of 
the MEA model is shown in Figure 2, and the structural element breakdown, 
numbering scneme, and connectivity are shown in Figure 3. 

The general SEA response equations for the six-element system with external 
excitation are: 
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N-f = number of modes resonant in element i 
Hi =* element i loss factor 

= power transfer coefficient for coupling between modes in 
elements i and j 

CO = center frequency of bandwidth 

Ei = mi total energy of element i 

( 1 ) 

Si = external acoustic or mechanical excitation in the bandwidth 
of interest 
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Figure 3. SEA Model Elements for Acoustic Test Configuration 




Since all of the elements are not physically connected to each other, the 
coefficient for power transfer between elements can be set to zero. 
This reduces the system of equations to: 
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2.2 DAMPING 

A loss factor of n = 0.01 vjas used for all other elements, based on Q = 100 
being a reasonable value for aluminum plates. The loss factor includes the 
contributions of internal energy losses, joint dissipations not included in 
the coupling loss factor, and losses due to radiation. 


2.3 MODAL DENSITY 

The number of modes (Ni) used in the response equations is determined by: 


Ni = ni(Af) 


where n^ = modal density of element i 
Af = bandwidth of the analysis. 

Element 1 of the MEA structure consists of the package support structure, 
thermal panels, radiator, interface support structure (ISS), and the signal 
distribution panel and components. The element modal density was calcula- 
ted by summing the modal density contributions of the element substructure. 
Since the thermal panels consist of fiberglass plates separated by insula- 
ting material, there is assumed to be no mechanical coupling between the 
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inner and outer plates (Reference 1). Therefore, they were modeled as two 
equivalent aluminum plates acting independently. The thermal panels were 
redefined as equivalent aluminum panels based on the respective material 
stiffness properties to aid in the ease of computations. All other parts 
are plates or are formed by multiple plate sections. Therefore, using the 
approximate equation for the high frequency modal density of plates, the 
element modal density could be found as: 

^ 2kp q 


where A 
t 
w 
E 
G 

V 

ni(f) 

Elements 2, 3, 4 and 5 consist of panels that are loaded by various compo- 
nents. Sevy and Earls (Reference 2) indicate that loaded panels will 
exhibit a greater stiffness and lower modal density than an identical 
unloaded panel. An increase in stiffness by a factor of 2 gave reasonable 
resi’l '.s in work done by Davis (Reference 3), and was assumed to apply for 
these elements. The resulting modal density for these elements could be 
found by the relationship 


^ 12w (1 


^ t 


) 


= surface area of plate 
= thickness of plate 

= weight density of plate (0.101 Ib/in^ for aluminum) 
= modulus of elasticity 
= 386 in./s^ 

= Poisson's ratio 
=6.30 modes/Hz 
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12w (1 -v=^) 


±_ 
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A 

t 


n 2 (f) = .0110 mode/Hz 
na(f) » .0076 mode/Hz 
04 (f) = .0237 mode/Hz 
ns(f) = .0083 mode/Hz 


Element 6 consists of the experiment mounting plates and the experiment 
assemblies. The modal density of the experiment mounting plates can be 
determined as for element 1. The experiment assembly is made up 0 ^ a 
cylinder capped by a dome. The approximate modal density of a cylinder 

is found by: 


ne(f) * 




for -^ < 1 
tr 


Substituting for kp(Cg,) as in the equation for element 1, the modal density 
for element 6 is expressed as 


n6(f) = 


2 / 2g . 
V 12w{l 


0/ 

7 12 w(l -v7 


la. 




■r 


< 1 


where fr = ring frequency of a cylinder (the frequency^whose wavelength 
equals the cylinder circumference) = J q 


The equation for the modal density of a doubly curved surface (Reference 4) 
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is an expression that is not readily evaluated, so the modal density of 
the experiment housing cap was approximated using the relation for a flat 
plate: 


n6(f) = Opiate ^cylinder + Ocap fo*" each experiment assembly. 
n6(f) = 0.609 mode/Hz 

= [0.444 (if . 

where fr = 3660 Hz. 

2.4 STRUCTURAL COUPLING 

All of the joints between elements are essentially plates joined at right 
angles. The relation for modal coupling for this type of joint is: 

A.. = iiL-iak-x 
Nj ZTT^fAi ^ 


•] 


f>i 

*r 


0.164 modes/Hz 




< 1 


for: 


Plate i 

Plate j 


S = l.OT^wCfct 


L = joint length 




hn* 


Di << Di 
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where 0 = plate rigidity = 

h-j = thickness of plate i 
h * 

hi < is the condition for Di « Dj 
i , j denote the elements that aire coupled 


The element coupling coefficients calculated using uiis relationship are 
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where <j>ij * <l>ji 
2.5 ELEMENT ENERGY 

The energy of the model elements is represented by 
Ei = m-j 
= mi 


<Vi^> 


<ai > 
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where < > indicates averaging over time and space 

= mass of element i obtained from MEA weight status summary 
mi = 615.14/g 
m 2 = 28.0/g 
m 3 = 108.84/g 

= 26.54/g 
ms = 670. 0/g 
me = 567.28/g 

2.6 ACOUSTIC POWER INPUT 

The external acoustic field is assumed to be reverberant, therefore the 
power input term can be represented by 

S, = 2tt^Cq~ Ai <P^> a Ni fsurfacpl 
w^(Au) mi 

The surfaces of the MEA package directly excited by the acoustic field 
were assumed to be comprised of the thermal panels, the radiator, and 
the interface support structure. The modal densities of these parts 
were found using the following relationships: 

n(f) = - 

2t! 

y i2w (1 - 

m = w A^ t/g 
N = n(f) Af 

It can be shown that 


/Mj.) , 4f 


A_\ A 

A_'j 1 

^ Surface 2 1 / 

wE H 

^ 'panels 

’^'radiator '^^/issj 


The thermal panels are made up of insulating material sandwiched between 
fiberglass cover panels. When two plates are separated by an air space 
or by loose filler materials (which is the case here), no mechanical 
coupling between the plates can be assumed. However, there is acoustic 
coupling, so the inner and outer panels were assumed to act independently 
with a resultant doubling of the thermal panel contributory surface area. 
Equivalent aluminum panels were calculated based on the ratio of the modulus 
of elasticity of aluminum to the modulus of the panel material. 

Therefore, 


= 7.412 X 10® Af 

^ surface 

The acoustic level criteria used in qualification testing of the MEA 
package are shown in Figure 4 and listed in Table 1. The values for each 
1/3 octave band center frequency were used as the acoustic pressure input 
where 


<p2> * 10^P^''^‘^(8.41 X 10'^®) 

Radiation efficiency values were obtained from Figure 5, based on a co 
incidence frequency Fc = 10,400 Hz (for an equivalent aluminum panel 
thickness tgq « 0.0368 in.}. These values are also listed in Table 1. 
The input term is therefore 

- Ca^ 10^P^''^°(8.41 X 10"^®) a(7.412 x 10®Af) 

= 

4tTf^ 


where Co = speed of sound in air 

= 13,400 in./s at 15“C, 60®F 


2-11 


Table 1 


SOUND PRESSURE 

LEVELS AND 


FOR USE IN 

FREp. 

SPL (dB) 

31.5 

120.0 

40 

122.0 

50 

124.0 

63 

126.0 

80 

128.0 

100 

130.0 

125 

132.0 

160 

133.0 

200 

134.5 

250 

135.0 

320 

135.0 ■ 

400 

135.0 

500 

135.0 

630 

133.0 

800 

131.0 

1000 

127.0 

1250 

126.0 

1600 

125.0 

2000 

123.0 

2500 

121.0 

3150 

119.0 

4000 

117.0 

5000 

115.0 


RADIATION EFFICIENCY 
SEA 

RADIATION EFFICIENCY 


.0191 
.0191 
.0191 
.0191 
.01 91 
.0191 
.01 91 
.0191 
.0191 
.0191 
,0191 
.0200 
.0209 
.0224 
.0240 
.0251 
.0275 
.0302 
.0355 
.0417 
.0501 
.0708 
.1059 


band mumbers 
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Figure 4 - One-Third Octave Band Acoustic Specificati 





2.7 INITIAL RESPONSE SOLUTION 


The response predictions for the elements were determined in each 1/3 
octave bandwidth from 31.5 to 5000 Hz by solving the system of equations 
for <Ti^>. The acceleration spectral density levels were then found where 


PS0(f)i 


— 2 
<ai > 


Af 


The root-mean-squared accelerations in the 1/3 octave band could be cal- 
culated from the relationship 

<Ti^> 



2.8 COMPARISON OF INITIAL RESPONSE PREDICTIONS AND MEASURED ACOUSTIC TEST 
DATA 

Comparisons between the SEA vibration response predictions and acoustic 
test data provided by MSFC were made by comparing plots of acceleration 
PSD's calculated for each structural element against plots of acceleration 
PSD's for measurements normal to the surface of the respective element. 

The acoustic test data consisted of 10 Hz bandwidth PSD measurements. In 
order to obtain a more appropriate comparison with the proper smoothing 
at higher frequencies, 1/3 octave band PSD's were formally computed from 
the 10 Hz bandwidth acoustic test data. The comparisons between the SEA 
predictions and all applicable test measurements are shown in Figures 6-11. 
The locations and orientations of the measurement points are identified 
in Table 2 and Figure 12. 

The predicted response levels for element 1 show very good correlation 
with measurements 24X and 25Y (Figure 6), except at low frequencies. Test 
points 24X and 25Y both measure accelerations normal to the surface of 
the structure, which is consistent with the direction of SEA predictions. 
The overprediction of SEA at lower frequencies can be related to the 
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Measurements for Element 2 








Comparison of Initial SEA Predictions and Acoustic Test 
Measurements for Elen»nt 6 





Table 2 


MEA ACOUSTIC TEST 
ACCELEROMETER LOCATIONS 


Accelerometer No. 

Axis 

Location Sketch No 

IX 


X 

Bottom Cold Plate and 

1 

2Y 

3Z, 

? Triaxial 

Y 

Z 

Battery Box Input 


4X 


X 

Battery Box Response 

1 & 2 


1 

X 

Primary Structure, 

3 

6Y 

> Tn axial 

Y 

Center Post 


7Zj 

1 

Z 



s 

I 

X 

Upper Cold Plate, Data 

4 


^ Triaxial 

i 

Y 

Acquisition Box Input 


lozj 

1 

Z 


nz 


Z 

Upper Cold Plate, Coolant 
Pump Input 

4 

12X] 

f Triaxial 

X 

Distributor Input, Vertical 

5 

13Y 

Y 

Panel 


14Zj 

1 

Z 


15Y 


Y 

Distributor 

5 


1 

X 

EAC Input 

1 & 6 

17Y 

► Triaxial 

Y 

18Z) 

1 

Z 



19X) 

1 

• Triaxial 

X 

Solenoid Valve Equipment 

6 

20Y> 

Y 

Panel 


2 izj 


Z 



22Z 


Z 

EAC Cap Response 

6 

23Y 


Y 

Equipment Panel Accumulator 
Input 

6 

24X 


X 

Thermal Panel #10 

7 

25Y 


Y 

Thermal Panel #6 

7 
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representation of all plate subelements as one tenn. At lower frequencies 
the individual plates have much lower modal densities than represented by 
the gross element model, which is a single "smeared" plate. 

For element 2, the SEA predictions are higher than the acoustic measurements, 
but the predicted response curve is seen to be converging to the measured 
value at higher frequencies (Figure 7). Measurements 12X and 14Z are 
parallel to the structure surface and are not appropriate for comparison 
with the SEA predictions. 13Y and 15Y measured accelerations perpendicular 
to the surface at two locations. When averaged they create a slightly 
higher resultant PSD. The slight difference in the prediction shape anr' 
the relative trends at higher frequencies seem to indicate that both 
frequency-dependent damping and coupling factor changes are required. 

The SEA prediction for element 3 is also higher than the acoustic test 
data. The general trends of the test data and the predictions are about 
the same except for a slight difference at higher frequencies (Figure 8). 

As in element 2, frequency-dependent damping and changes in the coupling 
factor should resolve the disagreement. 

The SEA prediction for element 4 matches up well with the acoustic test 
data (Figure 9). Again the prediction starts to fall below the data at 
higher frequencies, indicating variable damping is required. However, 
since element 4 is coupled through element 3, the indicated reduction in 
the predicted response for element 3 would cause a drop in the prediction 
for element 4. This seems to indicate that there is an alternate energy 
path not accounted for in the modeling. The alternate energy path could 
be due to additional acoustic excitation by the internal MEA acoustic 
field or through mechanical coupling to element 1 (photographs of the 
test article seem to indicate the latter possibility). 

In element 5, the SEA prediction is much higher than measurements 2Y, 3Z, 
and 4X located on the battery cold plate and battery box, and only slightly 
higher than data point IX (Figure 10). Measurements IX and 2Y monitored 
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motions in the plane of the structure and represent cross axis responses 
which do not apply to the SEA predictions. The convergent trend in the 
higher frequencies of the prediction and data from 3Z (normal to the 
battery mounting plate) indicates frequency-dependent damping is required. 
The generally high predictions indicate the coupling assumed in the analy- 
sis was excessive. 

The prediction for element 6 is much higher than the test data (Figure 11). 
Variable damping appears to be required because of the steeper slope of 
the prediction when compared to measurement 18Z (normal to the experiment 
plate). 16X and 17Y represent cross axis PSD's relative to the SEA pre- 
diction. Measurement 22Z appears to have a gain problem. Comparing data 
from 18Z and 22Z, the measured responses at lower frequencies follow the 
same trends except for a difference of two decades on the PSD (Figure 13). 
At lower frequencies the two responses should be nearly identical, because 
the stiffness of the canister structure dictates primarily rigid body mo- 
tions between the two measurements. At higher frequencies the shell 
response of the canister creates an additional departure from the SEA 
because it is not represented as an element in the model. 

The extra energy predicted by the SEA model for element 6 indicates that 
the coupling between elements 1 and 6 is not as large as predicted. The 
higher SEA prediction could also indicate that modeling of the experiment 
canister as part of the element was not correct. Treatment of the experi- 
ment and the experiment canister as mass loadings that reduce the modal 
density of the experiment panel could be more appropriate. This was the 
approach used in modeling elements 2, 3, 4 and 5. SEA response predic- 
tions may normally be expected to have good accuracy when there are more 
than 20 modes per analysis bandwidth (or above the ring frequency for a 
curved structure). Only element 1 meets this criterion in the 200-2000 Hz 
frequency range. 
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Section 3 


REFINEMENTS TO SEA MODEL 

The main refinements indicated by the comparison of the SEA predicted response 
levels with acoustic measurements were re-evaluation of the element damping 
and element-to-element coupling. Modifications to the modeling of the modal 
density of element 6 and identification of a possible alternative energy path 
to element 4 were also of interest. 

The modal density of element 6 was recalculated by assuming the experiments 
and experiment housings act primarily as mass loadings that stiffen the experi- 
ment mounting plates. Stiffening of the plates results in a lowering of the 
modal density of element 6. An increase in stiffness by a factor of 2 and a 
corresponding reduction in the modal density by 1//F was used in the calcu- 
lations for elements 2-5. The recalculated modal density of element 6 is 

ng = 0.080 modes/H2 

In the process of reviewing the structural drawings, one joint connecting 
element 1 to element 4 was found. This new path is diagrammed in Figure 14. 

The element coupling coefficient was calculated as explained in the section 
on structural coupling. 

i>i4 = 0.0008 

The system of equations describing the SEA model for acoustic excitation now 
becomes : 


ail 

ai2 

0tl3 

0tl4 

ais 

ai6 


El 


Si' 

®21 

tt 2 2 






E2 


0 

asi 


0t33 

CX34 




E3 

r = ^ 

0 

/ 

ai( 1 


0^4 3 

CK 44 




Ev 


0 

asi 




as 5 



Es 


0 

ag 1 





age 


He 

4 


0 . 


3-1 


< 3 * 



3-2 


Figure 14 - SEA Model Elements for Acoustic Test Configuration 



The need for frequency-dependent damping in elements 2-6 was indicated by 
the difference in slope of the predicted response levels and the test measure 
ments at higher frequencies. It was decided to decrease the rate of damping 
as a function of frequency above 250 Hz. The starting frequency of the 
decrement and the rate of decreasing damping were based on measurements from 
an acoustic test of Saturn IVB/V interstage panels (Ref. 3). The equation 
used to calculate the decreasing damping was of the form: 


where : 


f = frequency (Hz) 
n = element loss factor 

A graph of the damping values used in elements 2-6 is shown in Fig. 15. 

The basis for revising the element coupling coefficients was the observation 
that the outer element (element 1) prediction agreed with the test data, but 
the internal elements* response levels were overpredicted. The element 
coupling coefficients used in the analysis assumed the bolted connections 
between elements would be comparable to rigid joints. It is apparent this 
assumption was not valid. Methods of calculating the energy transfer through 
bolted connections were not found in any of the literature on SEA. It was 
therefore decided to arbitrarily reduce the coupling coefficient calculated 
for a rigid connection by an order of magnitude (10) for the bolted connec- 
tions. The energy path between the main resonator of element 1 (the thermal 
panels) and elements 2, 3. 5 and 6 passes through an intermediate subelement 
(the internal support structure-ISS). The thermal panels and elements 2. 3, 

5 and 6 are connected to the ISS by bolted connections. The "double" bolted 
joints in the energy transmission path seem to result in a reduction of the 
rigid joint coupling coeffieient by a factor of 100. 

The reduced coupling coefficients due to the "double" bolted joints are as 
follows; 
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<j>i2 = 0.0048//7 
<1.13 = 0.0116//7 
= 0.0008//7 
< 1.15 = 0.0093//7 
<1.16 ’ 0.020//7 
<i>34 ~ 38.6/v'^ 

The SEA predictions were recalculated using the frequency-varying damping 
and the reduced element coupling coefficients. 

3.1 COMPARISON OF REFINED SEA RESPONSE PREDICTIONS WITH MEASURED ACOUSTIC 
DATA 

The revised prediction response levels are compared with all applicable test 
measurements in Figures 16-21. 

The correlation between the predicted response levels and the test measure- 
ments remained very good for element 1 (Figure 16). The change in the system 
equations due to the revised damping and coupling was not enough to affect 
the prediction of levels in element 1 with any significance. 

For element 2, the SEA prediction now shows good correlation with the appli- 
cable test measurements (13Y and 15Y), except at lower frequencies (Figure 17). 


The revised SEA predictions for element 3 envelopes the measured data well 
and shows reasonably good correlation with measurements lOZ and IIZ (Figure 
18). The predicted response levels for element 4 with the addition of the 
coupling between elements 1 and 4 is shown in Figure 19. Since element 4 is 
also coupled to element 3, the reduction in the calculated response of element 
3 also lowered the response of element 4. The agreement between the measured 
data and the predicted response in the initial response calculations appears 
to be coincidental. The present predicted response levels are uniformly 
lower than the test data; however, the general trends of the SEA prediction 
and the test measurements match up very well for frequencies above 100 Hz. 
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This S66ms to lodicsts thot the use of frequency-dependent demping was correct, 
but that there is still an energy path not accounted for in the modeling. 

This energy path could be a mechanical path through the various hard lines 
connected to element 4 or through an acoustic path between element 1 and 
Element 4. 

In element 5, the SEA prediction matches 31 fairly well (Figure 20). The 
peaks and notches indicated in measurement 31 could be due to subelements 
that are modeled as mass loads in the SEA model. The largest test measurement 
(IX) is parallel to the surface of the battery plate and is not applicable to 
SEA predictions. 

The SEA response prediction for element 6 is still much higher than the test 
data (Figure 21). The use of variable damping adjusted the slope of the pre- 
diction so that the shape now matches that of the test measurements. The 
reduced coupling and reduced modal density has lowered the prediction levels, 
but more energy is still being represented in the model than is indicated by 
the test measurement. This could be explained by the location of the acceler- 
ometers in a corner of the experiment mounting plate. A more appropriate data 
comparison point would be one located in the center of the experiment mounting 
plate, which would be expected to register higher response levels. The data 
from measurement 22Z was dropped from consideration because of the suspected 
gain problem. 

In summary, the use of frequency-dependent damping seems to have been justi- 
fied. The trends of all the predictions match the test data at frequencies 
above 100 Hz. The reduction in coupling also seems to be justified because 
the predicted response levels match the test measurements except as noted in 
the above comparisons. 
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Comparison of Revised SEA Model Prediction and All 
Applicable Test Data for Element 5 









Section 4 


VIBRATION PREDICTIONS FOR PAYLOAD/SUPPORT STRUCTURE 
WITH COMBINED ACOUSTIC/MECHANICAL EXCITATION 


The main structure to be modeled is the MEA package modeled in Phase I 
of this contract (Vibration Predictions for a Pay load/ Support Structure 
with Acoustic Excitation). The Mission Peculiar Equipment (MPE) Support 
Structure has been attached to the MEA package for the SEA model for 
Phase II. Since the system of equations describing the model is linear, 
the predicted responses of the SEA model elements to mechanical excitation 
alone can be computed and added to the predicted responses for acoustic 
input alone to obtain a prediction for the combined acoustic and mechanical 
input response levels. 

The response of the MPE Support Structure to flight conditions was provided 
by MSFC. This known response was used as mechanical input to the system 
by the support structure. 


4.1 MODELING 

As in Phase I, the MEA package was divided into six elements. The MPE 
support structure was lumped into a seventh element. The structural 
element breakdown, element identification numbering, and connectivity 
are shown in Figure 22. 

For the seven-element system, the SEA response equations become: 
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Figure 22. SEA Model Elements for Flight Configurati 






ail 

0 ti 2 

Ol 3 

ai 4 

otl S 

ai 6 ai 7 


r 


fS:^ 

a 2 i 

O22 






E2 


S2 

as: 


O33 

a 34 




E3 


Sa 

a 4 i 


Oi,3 

a 4 4 



j 

> 

E4 

>=< 

S4 

asi 




ass 



Es 


Ss 

ae 1 





aee 


Es 


Se 

a? 1 





a ?7 






There is no external acoustic or mechanical excitation. 
Sx = S 2 = S 3 = $4 = Ss = $6 = $7 = 0 
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4.2 DAfviPING 

A loss factor of n * 0.01 was assigned to element 1. For elements 2-7, 
n was uniform below 250 Hz and frequency dependent above 250 Hz. The 
loss factor used in the prediction was calculated using the following: 

ri 2_7 ^ ^ < 250 Hz 

= 4.7687 * 250 Hz < f 
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4.3 MODAL DENSITY 


The modal density of elements 1 through 6 are the same as that used in 
the refined SEA model for acoustic excitation. 

ni = 6.84 modes/Hz 
nz - 0.0110 mode/Hz 
n 3 = 0.0076 mode/Hz 
n^^ = 0.0237 mode/Hz 
ns = 0.0083 mode/ Hz 
ns = 0.080 mode/Hz 

Element 7 basically consists of multiple plate sections. Therefore, the 
element density was found using the approximate equation for the high 
frequency modal density of plates. 

n? =• I t 

zi Ea 

y 12w{l -v2) 

= 2i683 modes/Hz 
4.4 STRUCTURAL COUPLING 

The element coupling coefficients calculated in refining the SEA model 
for acoustic excitation alone were used to describe coupling between 
elements 1 through 6. For the coupling between element 1 and element 7, 
reduction by an order of 10 because of the bolted joint was not used. 

The spacing of the bolts v«s not small enough to assume a connection 
approaching that of a rigid joint. Rather, a joint length of 3 inches 
was assumed around each bolt to be used in calculating the coupling co- 
efficient. The element coupling coefficients used for the mechanical 
input model are as follows: 
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<J>i2 = 0.0048//F 

«J>i3 = 0.0116//f 
<Pi>, = 0.0008//f 
(|)i5 = 0.0093//F 
<j>i6 = 0.020 //p 
^17 = 0.0044//P 
<))34 = 38.6 /y/f 

where 

4.5 ELEMENT ENERGY 

The element energy was handled as in Phase I. 

Where E^ = mi — ]-l 

mi =* 1.5936 lb-sVln. 

ni 2 = 0.0725 Ib-sVin. 

m 3 = 0.2820 Ib-sVin. 
niit “ 0.0688 Ib-s^/in. 
ms = 1.7358 Ib-s^/in. 
ms = 1.4696 Ib-s^/in. 
m? = 2.2021 Ib-sVin. 

4.6 RESPONSE SOLUTION 

The mechanical excitation input was provided by MSFC as X-axis, Y-axis, 

and Z-axis PDS levels (Table 3, Figure 23). The sum of the PSD levels’ 

was used as input in the SEA predictions. This was chosen since it is 
a more conservative approach than using the root of the sum of the squares. 
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Table 3 


MECHANICAL INPUT TO SEA MODEL 
Acceleration PSD (g^/Hz) 


FPvEQ 

X-AXIS 

(XIO'^} 

Y-AXIS 

(XIO’^) 

Z-AXIS 
(XI 0“^) 

X + Y + Z 
(XIO"^) 

(X^ + Y^ + 2^) 
(XIO"^) 

31.5 

o.n 

0.9 

0.15 

1.16 

.92 

40 

0.18 

1.45 

0.24 

1.87 

1.48 

50 

0.28 

2.25 

0.38 

2.91 

2.30 

63 

0.43 

3.6 

0.60 

4.63 

3.67 

80 

0.70 

5.8 

0.95 

7,45 

5.92 

100 

1.10 

9.0 

1.5 

11.6 

9.19 

125 

1.10 

9.0 

1.5 

11.6 

9.19 

160 

1.10 

9.0 

• 1.5 

11.6 

9.19 

200 

1.10 

9.0 

1.5 

11.6 

9.19 

250 

1.10 

9.0 

1.5 

11.6 

9.19 

315 

0.70 

9.0 

0.95 

10.65 

9.08 

400 

0.43 

9.0 

0.60 

10.03 

9.03 

500 

0.28 

5.8 

0.38 

6.46 

5.82 

630 

0.18 

3.6 

0,24 

4.02 

3.61 

800 

0.11 

2.25 

0.15 

2.51 

2.26 

1000 

0.07 

1.45 

0.10 

1.62 

1.46 

1250 

0.045 

0.95 

0.065 

1.06 

0.95 

1600 

0,029 

0.60 

0.040 

0.67 

0.60 

2000 

0.018 

0.36 

0.025 

0.40 

0.36 
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The response predictions for the elements 
bandwidth from 31.5 to 2000 Hz by solving 
<at^>. The acceleration spectral density 





were determined in each 1/3 octave 
the reduced set of equations for 
levels were then found where 


PSD(f)i = 


_ <ai' 


FAf 


RESPONSES^^ acoustic, MECHANICAL, AND ACOUSTIC/MECHANICAL EXCITATION 

The predicted MEA element PSD response levels for acoustic excitation, 
mechanical excitation, and combined mechanical /acoustic input are shown 
in Figures 24-30. They show the mechanical input has no effect on the 
SEA predicted response levels at lower frequencies, but becomes increasingly 
significant as the frequency gets higher. 



Figure 24 - SEA Response Predictions for Element 1 With Acoustic, 
Mechanical, and Acoustic/Mechanical Input 



Figure 25 - SEA Response Predictions For Element 2 With Acousti 
Mechanical , and Acoustic/Mechanical Input 







Mechanical, and Acoustic/Mechanical Input 





JjtA Response Predictions For Element 4 With Acoustic, 
Mechanical , and Acoustic/Mechanical Input 
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Figure 28 - SEA Response Predictions For Element 5 With Acoust 
Mechanical, and Acoustic/Mechanical Input 
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Figure 29 - SEA Response Predictions For Element 6 With Acoustic 
Mechanical, and Acoustic/Mechanical Input 
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Figure 30 - Acceleration PSD Level For Element 7 Used As 




Section 5 
CONCLUSIONS 

In the course of this phase of the study, conclusions were reached about the 
materials experiment assembly which was analyzed and about the application 
of SEA to payload vibration predictions. The following statements may be 
made, with regard to this prediction effort, which reflect observations on 
both points. 

1) The SEA predicted the response of primary structure to acoustic excitation 
quite well over a fairly wide frequency range when compared with the 
experiment. 

2) The SEA overpredicted the response of most internal structures when coup- 
ling factors were based on the existing body of SEA data. 

3) The SEA predicts that the contribution of mechanically induced random 
vibration to the total MEA response will not be significant except at 
higher frequencies. 

4) In most cases, two factors appear to be responsible for the overprediction 

a) Modal coupling factors at the joints were substantially lower than 
in the literature. 

b) Damping was frequency dependent. 

5) The power flow paths to the solenoid panel were not sufficiently definable 
to permit adequate resolution of the prediction errors. 

From these statements the general conclusion is that the SEA method has 
produced a prediction adequate for component specification purposes in all 
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but one case. Further work is needed to enlarge the body of data on internal 
losses, joint losses, and on modal coupling factors for various joint types, 
if accuracy is to be improved. In spite of this difficulty, SEA has the 
potential for replacing specific extrapolation methods for high frequency 
random vibration environmental predictions. 
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Section 6 

APPLICATIONS GUIDELINES 


1} Define response information required. 

2) Partition system into elements compatible with responses of Interest. 

3) Identify energy sources and energy paths between elements. 

4) Calculate 

i) Element damping values, 

ii) Element modal densities. 

iii) Element-to-element structural coupling coefficient, 

iv) Define acoustic and mechanical energy sources. 

5) Solve SEA system of equations for element energies. 

6.1 MODELING 

In applying SEA techniques, care should be taken in dividing the system 
structures into suitable elements. The number of elements and the sub- 
structure boundaries represented by the elements should be selected to best 
provide the desired response information. When choosing the model elements, 
two main limitations on the SEA technique should be considered. 

First, there is a tradeoff between the size of the structure represented by 
the element and the accuracy of the SEA prediction. This tradeoff is a result 
of the general proportional relationship between structure size and mode 
count. If there are many modes of an element excited in a frequency band, 
the energy of the system can be assumed to be distributed throughout the 
element and the modes within the analysis band. Averaging over the element 
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when this occurs then gives a valid approximation of actual response be- 
havior. Therefore the larger the substructure being modeled, the higher 
the mode count, and consequently the SEA element's average representation 
of the behavior of the substructure is improved. 

In order to obtain the high number of modes in each element, model elements 
are chosen as relatively gross portions of the structure. The finer details 
are “smeared” over the substructure in the SEA element model. Parts of the 
structure where response predictions are not required are lumped into other 
pieces of the structure for modeling. The averaging quality of SEA makes 
the lumping of structural parts possible since averaging over multiple parts 
of the structure is as valid as averaging over one part of the structure. 

As an added advantage, the lumping of substructures into one SEA element 
model results in a reduction of the bookkeeping required. 

The second limitation to be considered in modeling is the comparability 
of SEA predicted responses to the information required at discrete points 
on the structure. Since SEA averages over the entire element, the SEA cal- 
culated response is more comparable to' the behavior of some points of the 
structure, while the prediction is unrepresentative of the behavior of 
other points of the structure. For example, if a plate element is of 
interest, the response of the center of the plate would be roughly com- 
parable to the result of an SEA analysis while the response at a corner of 
the plate would not be comparable. The use of SEA and the interpretation 
of tne analysis results should be weighted by the two aforementioned 
considerations. 

The partitioning of the structure to determine SEA element boundaries and 
boundary conditions is made along actual structural joints or discontinuities. 
These divisions are made in accordance with the principles of SEA stated 
in the introduction. 

Once the structure has been partitioned, all energy sources should be identi- 
fied. The energy input to each element is usually identifiable as acoustic 
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or mechanical input. The definition of energy sources is used not only to 
indicate all sources of excitation of the system, but is also useful as an 
aid in limiting the size of the model. Any part of the structure that has 
a net energy flow of zero can be used as a model boundary. 

For example, consider the case where one element of the model is the external 
skin of an airframe. It is not necessary to include the entire external 
skin in the SEA model. The amount of external skin area to include in the 
model can be defined so there is a zero net flow of energy across the 
boundaries of this element of the model. Selection of the correct skin area 
yields a balanced system, with the energy flowing into the model subsystem 
mechanically from the remainder of the structure equal to the energy trans- 
ported mechanically out of the subsystem. In many cases this ideal condition 
can be approximately achieved by establishing the model boundaries at points 
halfway between major structural loading points, i.e., halfway between 
attach points of two adjacent equipment panels, halfway between the panel 
attach point and a fuel tank bulkhead, halfway between the panel attach 
point and a large component, etc. The effect on response predictions of 
incorrect estimation of the skin area -will be in essentially direct propor- 
tion to the error: selection of an area too large by 10% will result in 

predicted levels (g^/Hz) that are too high by approximately 10%. 

6.2 DAMPING 

A damping or internal loss factor for each element is required in the analysis. 
Although much work has been done in investigating damping, the choice of loss 
factor for each element is still fairly arbitrary. This arbitrariness is 
due to the loss factor representing not only the damping of the element but 
also the joint dissipations not included in the coupling loss factor and 
losses due to radiation. The best means of deriving the representation of 
the internal losses in the elements is engineering judgement based on exper- 
ience gained from testing similar structures. 

6.3 MODAL DENSITY 

The number of modes per frequency of an element is known as the element 


6-3 


modal density. This parameter is used to calculate the number of modes 
present in each SEA analysis band. Equations for calculating the approximate 
modal density of various structure shapes are available (Table 4, Ref. 4). 

Using these equations, the approximate modal density of each SEA element can 
be found. Since the structure represented by each SEA element is generally 
large, the element may actually be composed of hundreds of individual parts. 

The modeling cf all these parts by one SEA element is made possible by the 
averaging assumptions of SEA. A simple algorithm is used to calculate the 
SEA modal density parameter. The gross structure is decomposed into its 
various individual parts. These parts in turn can be decomposed, if necessary, 
into simpler components (i.e. an open box can be decomposed into five plate 
components). The modal density of these simple substructures can be calcu- 
lated using the approximate modal density equations. The modal density of 
the SEA element is then represented by summation of the modal densities of 
the various parts. 

One alternative to using the approximate equations is the packaged computer 
programs available that calculate the closed form solution of the mode 
count of various common structures (i.«. liquid-filled cylinders, pinned- 
end cylinders). A second alternative is the use of the results of low- 
frequency modal analysis tests. This can only be used where a low-frequency 
modal analysis of the structural element has been previously conducted. The 
results of the modal analysis are used to plot modes vs. mode number. If 
the trends of the graph are assumed to be valid at high frequencies, then the 
extrapolated slope is the value of the modal density of the element at the 
higher frequencies. 

6.4 STRUCTURAL COUPLING 

The structural coupling parameter is unique to SEA. The definition of this 
parameter is one of more important steps in the SEA procedure. However, 
except for formulations describing the rigid connection of a few beam and 
plate situations, very little information is available on methods of 
calculating this parameter (or on this parameter in general). This lack of 
knowledge is mainly due to the relative newness of the SEA method. The 
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Table 4 

MODAL DENSITIES OF SOME UNIFORM SYSTEMS 
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'(‘See next page for definitions of symbols. 
**N = Aw • n (w) 











Symbol Definitions for Table 4 



cross-section area 
surface area 

acoustic wave velocity 

Longitudinal wave velocity 

membrane wave velocity 

string wave velocity 

torsional wave velocity 

plate rigidity 
Young's modulus 
shear modulus 
thickness 

centroidal moment of inertia of A 
polar moment of inertia of A 
torsional constant of A 
length 

membrane tension force /unit edge length 

string tension force 

volume 

radius of gyration of A 

radius of gyration of plate cross section 

Poisson's ratio 
frequency (radians/time) 
material density 
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Mu»t1ons tor calculating the coupling of the various beam and plate Joints 
a« found in References 1, 6 and 7. Coupling relations for a tew Wtcal 
joints found in testing are also available in the literature. 

When structural Joints or energy paths (i.e. hard 

„he any previously defined in the literature are ^ - 

the structural coupling must be empirically assigned to that joint. The 
Ln coupling parameters can sosmtinms be used as gu^elines oc u^per . 
aolutions, but the reduction required for the approximate 
the structure is mainly a matter of Judgement based on experience. 

6.5 ACOUSTIC INPUT 

The acoustic coupling representation used assumes that the external acoustic 
field is reverberant. The expression is of the form: 

Cj A^ <p2> a 
S 

(Auj) jn. 

If other acoustic fields are to be used as input, an "equivalent reverberant 
field would need to be defined. 

The radiation efficiency term. o. may be approximated using procedures out- 
lined in Reference 4, or may be determined from Figure 5. 

6.6 MECHANICAL INPUT 

A representation of mechanical input to the SEA model can be used when 1) 
thl resp^nse of one element of the system is known, or E) when the system is 

coupled to a known source of excitation. 

In case 1 the known response of the element can be plugged directly 
! stem of’equations. This reduces the order of the problem by one because 
it reduces the number of unknowns in the system of equations to e so 

by one. 
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I , In case 2, the known external mechanical energy source is incorporated into 

the SEA model as another element. The system of equations is then handled 
■| as in case 1. Since the energy level of the courcc is known, the size of 

R the problem is not affected by the addition of the mechanical source as 

another element. The systan of equations can then be solved for the other 
element energy levels. 

6.7 SEA RESPONSE PREDICTION EQUATIONS 

The SEA response prediction equations consist of a set of equations that 
describe the energy state of element of the model. Each equation describes 

1) the power into the element from an external source 

2} the power dissipated within the element 

3) the net power transferred between the element and any other 

element or elements coupled to it. 

The general system of equations is of the form: 



N^ = number of modes resonant in element i 
= element i loss factor 
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* power transfer coefficient for coupling 
between modes in elements i and j 
(symmetric 

o» = center frequency of the bandwidth 

E,- = = total energy of element i 

1 * 2 
ii)'~ 

S. = external acoustic or mechanical excitation 
in the bandwidth of interest 

The unknowns of the set of equations are the total energy at the average 
frequency m. The system of equations are solved simultaneously for the un- 
known energy levels. 

A two-element system can be used to see how these equations are divided. 



'a 

n 

Sb 

1 

_:.a. b ^ 

1 

cr 

1 : 

L_i 


°a D, 

' 

b 


S = power introduced into element a from an external source in the 
bandwidth of interest. 

D = power dissipated within element a in the bandwidth of interest. 

cl 

b * power transmitted from element a to element fa (= -P|j^ 
in the bandwidth of interest. 

The energy passing through the two elements in a single bandwidth can be 
expressed as: 
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^ *^a,b ^a 

\ * ^b,a ^ h 


To calculate solutions over the spectrum of interest, the system of equations 
is solved for each contributing frequency band. 

The energy dissipated per unit time is defined in terms of the element loss 
factor as 

“a " “"a^a 

where 

o> = angular frequency (average) of system 
n. * element a loss factor 

a 

E, = total energy of element a 

a 

The net power transmitted from the resonant modes of element a to the resonant 
modes of element b is 

’’a.b ” *'*b'^a,b^a “ ^0*^0, b^b 

= (power transmitted from b to a) - (power transmitted 
from a to b) 


where 


N, * number of modes resonant in element a 
a 


(j»a,b = power transfer coefficient for coupling between modes 
through the structural joint 

Performing the indicated substitutions, 

“’’a^a * ^b^a,b^a " '*a*a,b^b * h 

“'’b^b ^a*a,b^b " \*a,b^a “ ^b 
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Regrouping the terms of the equation the general form of the response 
solution for two elements is obtained. 




-"l, *b. 


-"a ♦ab 


(“”b * “a *ba> 


n, N» <? and S are defined parameters, and the set of linear simultaneous 
equations can be solved for the unknowns Ea and Eb at the average frequency 

aj. 
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